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Abstract 



The exphcit solution is obtained for four-wave mixing UJ4 = uJi — uj2 + 
of two strong fields Ei,E^ and two weak fields E2,E4 in a four- level system 
with the large Doppler broadening. The resonance of the mixing coefficient 
dependence on intensity is found around Eidi = E^d^, where di^s are the 
dipole moments of corresponding transitions. The effect is interpreted as an 
intersection of quasi-energy levels. Up to 6 peaks appear in the dependence of 
conversion coefficient on the detuning of the probe field £'2- An unexpected 
additional pair of peaks is a consequence of averaging over velocities and 
disappears at low temperature. The results allow us to interpret saturation 
behavior in recent experiments on the mixing in sodium vapor. 
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Typeset using REVT^ 



Four- level system is a promising object for resonant optics and spectroscopy owing to a 
great variety of nonlinear effects. They are nonlinear interference, inversionless gain, res- 
onance refraction, electromagnetically induced transparency, optically induced energy-level 
mixing and shifting, population redistribution etc (see |jl|J^ and citations therein). Recent ex- 
periments on continuous four-wave frequency mixing of Raman type with sodium molecules 
in a heat pipe gave interesting dependencies of generated wave power on frequencies 
and intensities of the incident waves. In particular, the dependence of output power on the 
first strong field intensity was saturated in experiment in thin media, whenever that on 
the third wave intensity demonstrated the linear growth. The measurements were taken at 
large Doppler broadening while the nonperturbative theory is proposed for atoms at 
rest. The saturation of four- wave mixing efficiency as a function of intensity was observed 
in optically thick lead vapor interacted with megawatt pulses 0. In case the Rabi splitting 
exceeded the inhomogeneous width. 

Development of nonperturbative theory from the mathematical point of view involves the 
solution of the set of 16 algebraic equations for steady-state elements of atomic density matrix 
for four-level system. The problem is only of analyzing the resultant awkward expression 
and to average it with Maxwellian distribution over velocities. In present paper we study the 
particular case of alternate two strong and two weak fields interacting with 4-level system 
having some symmetry. The 4th degree equation can be reduced to biquadratic one, then 
the integration is possible analytically. 

A simple explicit formula for nonlinear susceptibility at zero frequency detuning displays 
the resonant behavior as a function of a strong field at fixed another strong field. We 
interpreted it as the coherent effect owing to the intersection of quasi-energy levels. The 
susceptibility as a function of weak field frequency a;2, has 6, 4 or 3 peaks. The profile also 
displays the important role of the Rabi splitting. 

Let us consider the conversion of two strong incident waves Ei^^ resonantly interacting 
with opposite transitions gl, mn and the weak field E2 near the resonance with transition 
gn into the 4-th output wave E^, inset of Fig. 1. The electric field in the cell is 

E{r, t) = ^E^ exp (iuJut - ik„fj , (1) 

where Ei, is the amplitude of i^-th field, cu^, ky are the frequency and wavevector. The strong 
fields are also near resonance, oji ~ ojguuj?^ ~ cUmn, where Uij = {Ei — Ej) /Ti are transition 
frequencies between energy levels Ei and Ej. The intensity of the 4-th wave, that appears 
during the process of mixing, is also being small. Its frequency and wavevector satisfy the 
phase matching condition 

uJi = uoi - UJ2 + t^3, ki = ki-k2 + ki. (2) 
Maxwell equation for the output wave can be reduced to 

d£^4 2'KiuJrnldml 



dx 



{Pml) , (3) 



where x is the coordinate, dmi is the matrix element of the dipole moment operator d, c is the 
speed of light, pmi is the coherence at transition ml, angular brackets denote the averaging 
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over velocity distribution. We should calculate pmi as a function of input amplitudes -E'1,2,3, 
their wavevectors /ci,2,3 and frequency detuning of the weak field. 

For this end we solve the equation for Wigner's atomic density matrix (see 0) 

where v is the atomic velocity, jij are relaxation constants, Qj = exp(— '5^/f|.)/f|.7r'^/^ 

is the Maxwellian excitation function, V = —E{r,t)d/2h is the operator of interaction, 
i,j = m,n,g,l. 

To the zeroth approximation we can neglect both the weak fields -E'2,4 0. The set 
boils down to finding out populations pj = pjj and coherences pi = pgiexp{—iuJit + ikir), 
P3 = exp(— ituat + i/c3r) of a pair of separated two-level systems. The solution is written 

as 

pi = Ni- Ni- 



pi 



P3 



^7„(r23 + (]f)' 



^g3iV3r^ 



where Nj = qj/'jjj,j = m,n,l,g are the unperturbed populations, Ni = Ni — Ng, N3 = 
Nn — Nm are the population differences at "strong" transitions, 71 = 'Jgi,'j3 = ■jmn are their 
homogeneous width, Gi = Eidgi/2h, G3 = E^dmn/'^ft- are the Rabi frequencies, fl^ = Q^—k^v 
is the Doppler-shifted detuning Qi = ui — Ugi, = — Umn, = lu + 'i^l, 

r^, = 7? + 2|Gi|V(7;^+7r'), 
r^3 = 7| + 2|G3r73(7^'+7„-') 

are the homogeneous widths including the power broadening. 

Weak fields with amplitudes G2 = E2dgn/2h, G^^ = E4^dmi/'2h lead to appearance of cross- 
coherence between levels belonging to the opposite two-level systems p2 = Pgn, P4 = Pmi at 
the allowed transitions, as well as at the forbidden transitions ps = Pgm, Pe = Pni- To the 
first order one can neglect the infiuence of these fields to the populations. The set of 4 
algebraic equations appears for the nondiagonal matrix elements: 
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T2P2 - iGipl + iGsp^ = -iG2{pg - Pn), 

Tlpl + iGlpl-iGlp5 = iGl{pm- pi), (7) 
Tsps - iGipl + iGlp2 = iG2p*3 - iGlpi, 
r*pg + iGspl - T'G*ip2 = -iG2p*i + 'i'Glps- 

Here 72 = 7gn,74 = Imi are the constants of relaxation of the coherence at the allowed 
transition, 75 = 7^^, 76 = Ini are the constants for forbidden transitions, Q'^ = Q[ — Q'^, Q'q = 
Q[ — Q2 are frequency detunings at these transitions. 

The solution of Eq. (^) for the nondiagonal element at output transition ml can be 
presented as 

p* = -if34GlG2Gl - ia^Gl (8) 

During the initial step of the mixing the generated field is small, IG4I <^ 1^21, that enables 
one to neglect the absorption and to find the coefficient j3/i only. We found intensity of 
output wave within the thin medium approximation by the integration of Eq. (0) from a; = 
to the length of cell L 

h{L) = '^1^^^ (A) {dgiei){dgne2) x 

{,drane'i){draiei) hhh, (9) 

where e^, is the polarization of z/-th wave, I,j is its intensity. We find coefficient (3^ comparing 
Eq. (|^ to solution of the form (||). 

/34 = {{T, + r*)(p, - p„) - §;{\G,? - IGar 

-r2r5)- J^(|G3r-|Gip-r2r*)^. (lo) 

Here elements pg, pn, pi, P3 are defined by Eq. (|^), (|^), the determinant of set (|^) is 



d = T2T,tit; + {\Gi\^-\G3\Y 

2' 



+^(|G'ir + |G'3|')(r2 + r:)(r5 + r 



-liiGii'-m'm-TDiT.-Ti), (11) 

the polynomial of 4th degree in velocity. The averaging of coefficient jS^ over velocity is 
possible by the residues theory for the Doppler limit vt ^ 00. 

To examine the intensity dependence of coefficient let us consider the case of equal 
relaxation constants of the levels jj = 7,J = g,l,m,n, excitation of the lower level only 
Qj/j = NSji, the resonant strong field detunings fli/ki = Q^/k^ vt, and equal wavenum- 
bers of both weak fields \k2 — k^] -C {k2kik5kQY^^ . The last condition is natural in down- 
conversion scheme. In view of phase matching condition (Q) it is felt that the weak field 
detunings depend on single parameter f2: 0.2 = k2^i/ki + Q,Q4 = k^S^i/ki — VL. Within the 
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assumptions one can also see that = k2,ke = k^. If all the wavevectors are parallel, then 
the expression for {j34) assumes a simple form 



dx 



\Gi\ 



^vt'^ J D{x)r^,^ + klx^' 

— OO 

C{x) = A\Gi\^iz + (7 - ikix) x 
31 - {1 - i{k2X - Vl){'-f - i{kr,x - ^l)) 



(12) 



Here x = k2v/k2 — ^li/ki, z = Q — i'j, = 7^+4|Gip is the saturated width. Determinant 
D(x) turns to be a function of 



D{xj = K^X 
.2 In i\^2 
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k2k^ 



(13) 



Ai = (//2-l)^^-|Gi| 



^2 



z^-{\G,\ + \G;\f 



z'-{\G^\-\G,\r 
The detuning dependence of IA2I takes the minimal values at 

VL = ±\Gi\ ± IG3I. 

It is a consequence of the level splitting by the strong driving field. Note that at |Gi 
two points of minimum merge together. The reason is equal Rabi splitting for each 
The simple form of the determinant ([T3|) allows calculating mixing coefficient ( 
plicitly 



(14) 

= ICsl 
level. 
IT2D ex- 



7 + iz^ 



+ 



4^216-112 + 7(;22 ^ IG3 



R 

R Tsi 



(15) 



where R = y 2(A2 — Ai), ^R > 0. The branch of two-valued function A2 should be chosen 
according to the following rules 

3?A2 < at P+ < \n\, 3f?A2 > at \n\ < P^, 
sign(^5A2) = sign Q at P_ < \Q\ < P+, 

where P± = ||Gi| ± IG3II. 

The mixing coefficient | {P4) p calculated from Eq. (|15|) is plotted in Fig. 1 ( 
function of detuning Q. The coefficient has 4 peaks at points given by ([T^) . At equal 
distances between quasi-energy levels |Gi| = {G^l two central peaks coalesce in the center 
f2 = 0, Fig.l (c). Except of the zeros of A2, zeros of R{fl) may add two peaks near the 
center, Fig. 1 (b). The additional central peaks are absent for motionless atoms since only 
four transitions are possible between two pairs of splitted quasi-energy sublevels. These 
peaks are contrast at Gi > 6*3 ^ 7 and disappear at IG1/G3I < /i/V/i^ — 1. 

The value | (/54) p at the exact resonance fij^ = 0, z/ = 1, . . . , 4 is shown in Fig. 2 as a 
function of |Gip. The sharp peak at |Gi| = IG3I confirms the qualitative interpretation 
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of the effect as tlie intersection of quasi-energy levels. Inset in Fig. 2 illustrates why the 
maximal conversion occurs when the Rabi splitting in opposite two-level systems are equal. 
Here the cross-transition from the upper sublevel of level m to the upper sublevel of level 
n has the same frequency as the transition between their lower sublevels. In this case only 
3 resonances remain in the spectrum, Fig. 1 (c), with the overpowering maximum in the 
center. The resonance condition = IG3I brings the maximum conversion efficiency in 
the intensity dependence. 

The splitting effect is evident from experimental results on resonant four wave mixing in 
Na2 . The main feature is the saturation of output power as a function of one strong field. 
The conditions of experiment are generally satisfy the above model: (1) down-conversion 
level scheme 074 < ui (see inset, Fig. 1) with kiVT = 7.0 • 10^ s~\ k2VT = 6.5 ■ 10^ s~^, k^vr = 
5.2 ■ 10^ s~^, k^VT = 5.7- 10^ s~^; (2) all incident waves 1,2,3 are generated by external lasers; 
(3) the region of interaction is short enough (nearly 1 cm), the model of thin media can be 
treated; (4) estimated level parameters are A*"; ~ lO"*^^ cm~^ ^ Nn ~ 10^^ cm~^ ^ A^^, A^^, 
7^ ~ 7g ~ 2 ■ 10® s~^, 7„ ~ 7; ~ 2 ■ 10^ s~^. Slightly noncoUinear geometry (mixing angle 
6 ~ 10~^) leads to an effective broadening Auj ~ kvT ■ ~ 10® s~^. Another factor is usual 
jitter of laser frequencies, especially for dimer and dye lasers, Auj ~ (2-h4)-10® s~^. Thus, the 
effective value 7 = (3 -i- 6) ■ 10® s~^ seems reasonable; (5) the maximal field values estimated 
from the focusing geometry jdr^^ ~ 10^ s-\ jGaT'''' ~ 2 • 10® s-^ IG'al'"^^ -5-10® s'^ 
nearly correspond to the condition of two strong fields. 

The resonance condition IGil = \Gz\ may result in peaks as in dependence /34(/i), as 
in dependence I3i{h)- If > IG'sl™'^'', the peak is seen only in [3i{Ii). The width of 

the peak is determined by the decay rate 7. Since in the experiment 7 ~ 1^3 1, the peak 
is wide. Fig. 2(b), and gives a smooth saturation curve /4(/i), Fig. 2(c). According to this 
consideration the saturation of in the experiment (boxes in Fig. 2) is observed at 

ICil > IG3I and there is no saturation for I^ll^). Note that such behavior was observed 
for different values of I2 varied by one order. Under the opposite experimental condition 
1^ < |G3|™'^^ the dependencies IaIIi) and I^lls) change their behavior in agreement 

with the consideration. 

Thus, the model explains quantitatively the main features of the measured saturation 
curves. To observe the sharp resonances arising from Rabi splitting the stabilization of laser 
frequencies seems to be important. To increase the efficiency of conversion into the 4th wave 
it is necessary to tune up the laser frequencies to corresponding peaks. The optimum at 
fl^ = corresponds to equal Rabi frequencies \Gi\ = IG3I. 

Authors are grateful to S.G. Rautian, A.M. Shalagin, and M.G. Stepanov for fruitful 
discussions, B. Wellegehausen and A. A. Apolonsky for clarifying the details of experiments. 
This work was partially supported by Deutsche Forschungsgemeinschaft, grant WE 872/18- 
1. 
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LIST OF CAPTIONS 

Fig. 1. Conversion coefficient | {j34) p (arb. units) as a function of detuning Q of the 
second field at [Gil = 1, IG3I = 0.5, = 7.0, fcafr = 6.9, 7 = 0.2 (a), 7 = 0.02 (b), and 
7 = 0.02 at |Gi| = IG3I = 0.5 (c) (all frequencies are in ns^^). Inset is the level diagram of 
four-level system interacting with two strong driving fields at the opposite transitions (solid 
arrows) and two weak fields (wavy arrows). Dotted lines show the forbidden transitions. 

Fig. 2. Conversion coefficient | {P4) p (arb. units) vs at {G^l = 0.5, = 0, 

k^VT = 7, k2VT = 6.5: 7 = 0.06 (a), 7 = 0.6 (b), and vs at 7 = 0.6 (c). The 

parameters for (b), (c) correspond to experiment, all frequencies are in ns^^. Boxes denote 
the experimental points from 0. The inset illustrates the Rabi splitting of dressed states. 
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